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ABSTRACT 

We present a long-term X-ray flux and spectral analysis for 1RXS J170849. 0—400910 using 
Swift /XRT spanning over 8 years from 2005-2013. We also analyze two observations from Chandra 
and XMM in the period from 2003-2004. In this 10-yr period, 1RXS J170849. 0—400910 displayed 
several rotational glitches. Previous studies have claimed variations in the X-ray emission associated 
with some of the glitches. From our analysis we find no evidence for significant X-ray flux varia- 
tions and evidence for only low-level spectral variations. We also present an updated timing solution 
for 1RXS J170849. 0—400910, from RXTE and Swift observations, which includes a previously unre- 
ported glitch at MJD 56019. We discuss the frequency and implications of radiatively quiet glitches 
in magnetars. 

Subject headings: pulsars: individual (1RXS J170849. 0—400910) — stars: neutron — X-rays: general 


1. INTRODUCTION 

Magnetars are a type of pulsar that exhibit exotic and 
often violent properties. Their defining characteristic is 
that they are powered not by their rotation, as are most 
Crab-like pulsars, but by the decay of their high mag- 
netic fields. Because of the energy provided by the mag- 
netic field decay, their X-ray luminosities are generally 
higher than their rotational spin-down energies. They 
often display outburst activity during which they can in- 
crease their brightness by an order of magnitude or more 
and emit short (~ 10ms to ~ Is in duration) energetic 
bursts. Previously, magnetars had been classified into 
two observational categories: Anomolous X-ray pulsars 
(AXPs) and Soft-gamma repeaters (SGRs). However, 
these two “classes” appear to be merely different ends 
of the magnetar beha vioral spectrum. For a review see 
IWoods fe Thompsonl (|2006f ). 

Glitches in magnetars have been observed both with 
and without associated radiative chan ges. Out_ of th e 26 
magn etars and magnetar candidate!! (lOlausen fc Kasnil 
12013 ) only five are monitored sufficiently frequently 
to detect unambig uously the occurance of glitches 
(|Dib fe Kaspil 120131) . Of those five, one, IE 1841—045, 
has never displayed any ra diative activity associ- 
ated with its glitches (|Zhu fc Kaspil 12010( 1 whereas 
IE 1048.1-5937, IE 2259+586, and 4U 0142+61 have 
had radiat i ve ev e nts during some o r all of their g l itches 
(IDib et al.l 120091: Kaspi et al.1 120031 iGavriil et al. 120111 : 
iDib fe Kaspill2013 ) . It is important to determine whether 
or not there is a generic connection between magnetar 
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glitches and radiative events because it can help us de- 
termine the physical orgin of these phenomena. It seems 
reasonable that magnetospheric mechanisms, because of 
their external nature, are likely to be accompanied by ra- 
diative changes whereas internal mechanisms could pro- 
duce radiatively quiet glitches. 

1RXS J170849. 0—400910 (hereafter referred to as 
RXS J1708 for brevity) was first ide ntified as an X-ray 
source in the ROSAT all-sky survey flVoges et al.ll 19991). 
It wa s first discovered as a pulsar by iSugizaki et al.1 
(fl997T) , using ASCA data, who suggested that it was an 
AX P based on its X-ray spectrum and 11-s spin period. 
Ilsrael et ~ah] (11999T) measured a period derivative typical 
of AXPs for RXS J1708, confirming that the source is an 
AXP, and thus a magnetar. 

RXS J170 8 was the first magnetar observed to glitch 
(|Kaspi et al.1 1200dh . It has since been found to glitch 
several more time s (Ilsrael et al.l 120071 : IDib et all 120081 : 
IDib fe Kasoil 1201 31) . Note t hat some of the glitches re- 
ported in Ilsrael et al.l (1200711 a re considered to be glitch 
candidates in IDib et al.l (12008H as they c ould be consis- 
tent with timing noise. IRea et al.l (|2005T ) first suggested 
that RXS J1708 exhibited post-glitch X-ray flux variabil- 
ity based on a 2003 XMM observation. They reported 
an XMM flux that was significantly lower than preceding 
Chandra and BeppoSAX observations. Futher evidence 
for flux variability was claimed ba sed on additio nal Swift 
and INTEGRAL observati ons JCampana et al.1 120011 : 
IGotz et ail 2 0071: Ilsrael et alJl2007li . However, puzzlinglv. 
variability at the level claimed in these studies was not 
seen the pulsed coun t rate as measured by frequent obser - 
vations with RXTE (IDib et al.ll2008t IDib fe Kaspill2013H . 

In this paper we analyze all the available Swift X-Ray 
Telescope (XRT) data from the period of the claimed 
variability to present day. We also use one XMM and 
one Chandra observation that were performed prior to 
the start of the Swift observations. We then use the 
measured spectral and flux values to constrain the level 
of source variability. We also present an up-to-date 
timing soluti on which continues the RXTE timing of 
IDib fe Kaspil (12013T) using Swift. We then discuss the 
occurance of radiatively quiet glitches in magnetars. 
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2. OBSERVATIONS 
2.1. Swift Observations 

RXS J1708 was observed by the Swift XRT fre- 
quently between 2005 and 2010. Beginning in July 2011 
RXS J1708 was observed as a conti nuation of the RXTE 
timing campaign summarized by iDib fc Kaspil (2013). 
Here we use all available archival Swift data in that 
time period in both Windowed Timing (WT) and Photon 
Counting (PC) modes. There were 80 observations for a 
total exposure time of 268 ks. Table [l] shows a summary 
of the Swift observations used in this work. 

We downloaded the unfiltered Level 1 data from the 
HEASARC data archive and ran the standard Swift data 
reduction script xrtpipeline using the source position 
of 17 h 08 m 46.87 s , -40° 08' 52.44" (llsrael et al.ll2003li and 
the best available spacecraft attitude file. Events were 
then reduced to the solar-system barycenter using the 
same position. For WT mode, a 30-pixel long strip cen- 
tered on the source was used to extract the source events 
and a 50-pixel long strip positioned away from the source 
was used to extract the background events. For PC mode 
observations, an annular region with inner radius 3 pix- 
els and outer radius 20 pixels was used. The inner region 
was excluded to avoid pileup of the source. An annulus 
with inner radius 40 pixels and outer radius 60 pixels was 
used as the background region. 

For WT mode data, exposure maps, spectra, and an- 
cillary response files were created for each individual or- 
bit. The spectra and ancillary response files were then 
summed to create a spectrum for each observation. For 
the PC mode data, exposure maps, spectra and ancillary 
response files were created on a per observation basis. We 
used response files for spectral fitting from the 20120209 
CALDB. 

The use of exposure maps when creating the ancillary 
response files is especially important for Swift data, as 
there are columns of bad pixels which can disrupt the 
PSF of the source for parts of certain observations. Or- 
bits were not used in the observation if the bad columns 
were found to be within 3 pixels of the source position. 

For WT mode data we selected only Grade 0 events for 
spectral fitting as higher Grade ev ents are more likely t o 
be caused by a background event dBurrows et al.ll2005ll . 
In PC mode we used the standard Grade 0-12 selection. 

2.2. Chandra and XMM Observations 

In this study, we also reprocessed archival data 
taken with the Chandra X-ray Observatory and XMM- 
Newton. To avoid pileup, we used the Chandra 
continuous-clocking (CC) mode observation (ObsID 
4605) and the XMM PN small-window mode data (Ob- 
sID 0148690101). The former was taken on 2004 July 3 
with the ACIS-S detector in CC mode, which has a time 
resolution of 3 ms. The total exposure was 29 ks. The 
XMM observations were made on 2003 August 28. The 
PN and MOS detectors were run in small and large win- 
dow modes, with 0.5-s and 6-ms time resolution, respec- 
tively. As the source is bright, the low time resolution of 
the MOS data results in significant pileup. Therefore, we 
focused only on the PN data. After filtering for periods 
of high background, we were left with 35 ks of exposure. 
This is equivalent to 24 ks of live time since the small- 
window mode has an efficiency of 70%. 


We processed the Chandra and XMM data using CIAO 
4.4 and SAS 11, respectively. The source spectrum was 
extracted using a 6"-wide region from the Chandra ob- 
servation and a 40"-radius aperture from the XMM PN 
data. For the Chandra observation, the background spec- 
trum was extracted from the entire ID CC-mode strip 
excluding the inner 1' closest to the source. For the 
XMM observation, the background spectrum was ex- 
tracted from two 40"-radius circular regions placed away 
from the source. 

3. ANALYSIS & RESULTS 
3.1. Flux and spectra 

We first fit the spectra for each individual observation 
with a photoelectrically absorbed power-law model. The 
spectra were fit with a single Ah usi ng t he XS PEC tbabs 
model with abundances from IWilrns et a hl (12000 ). and 
photoelectric cross-sections from lYerner et al.l ( 199611 . 
We used XSPEC0ittp://xspec. gfsc.nasa.gov with Cash 
statistics (lCashlil979l) to fit the spectra because of the 
low number of counts in the Swift observations. The 
grey points in Figure [l] show the results of the spectral 
fits to the individual observations. The typical uncer- 
tainties in the spectral parameters vary widely due to the 
large range in exposure times. In order to place the best 
constraints on the variability, we consider PC and WT 
modes separately. This is because the two modes are cal- 
ibrated to within only 10% of each other (A. Beardmore, 
private communication). The mean and standard devia- 
tion of the 1-10 keV absorbed flux for the PC mode data 
are 4.0 x 10” 11 erg cm -2 s -1 and 1.9 x 10” 12 erg cm -2 s _1 , 
respectively. The PC mode photon index has a mean and 
standard deviation of 3.1 and 0.08. For WT mode the 
mean and standard deviation are 4.0 x 10 -11 erg cm -2 s _1 
and 2.1 x 10 -12 erg cm -2 s -1 for flux and 3.2 and 0.07 for 
the photon index. 

In order to better constrain the variability, we then sep- 
arated the Swift spectra into sets of observations nearby 
in time (see Table [l]). Within each set, the 1-10 keV flux 
and photon-index were consistent with being constant 
(i.e. the y 2 values of fits to a mean value in each set were 
consistent with being drawn from a y 2 distribution) . We 
fit the sets of Swift observations as well as the XMM and 
Chandra spectra with a photoelectrically absorbed power 
law. Each Swift set was fitted with the same model with 
all spectral parameters the same from observation to ob- 
servation within the set. All parameters were allowed 
to vary from set to set except for Ah which was tied to 
the same parameter for all sets and was measured to be 
(2.434 ± 0.008) x 10 22 cm -2 . We did not fit the conven- 
tional but more complicated blackbody plus power-law 
model because the addition of the extra blackbody com- 
ponent did not improve the goodness-of-fit significantly 
for any of the Swift sets. Although additional compo- 
nents are significant for the XMM and Chandra spectra, 
we opted to use a single component model because only 
Swift d ata a re used here to constrain the variability (see 
Section E|) and using a single-component model sim- 
plifies the comparison of spectral properties. The joint 
power-law fit to the sets of observations provided a Cash 
statistic of 32806 and a Pearson y 2 of 37191 for 33571 
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degrees of freedom. This corresponds to a reduced y 2 of 
1 . 1 . 

Figure [T| shows the 1-10 keV absorbed flux and power- 
law index as a function of time resulting from the spectral 
fit to the sets. The mean and standard deviation of the 
flux for the PC mode sets are 4.0 x 10” 11 ergcm^ 2 s _1 
and 8.4 x 10~ 13 ergcm^ 2 s -1 , respectively. The PC mode 
photon index has a mean and standard deviation of 3.1 
and 0.07. Thus for PC mode, the maximum variability 
allowed by 3<r confidence intervals (3 times the standard 
deviation divided by the average value) is 6.3% for the 
flux and 7.5% for the photon index. In WT mode, the 
mean flux is measured to be 4.1 x 10 -11 ergcm~ 2 s -1 
with a standard deviation of 5.9 x 10 -13 ergcm^ 2 s _1 for 
a maximum variability of 4.3%. For the photon index, 
the mean and standard deviation are 3.2 and 0.04 for a 
maxiumum variability of 4.0%. Compared to the stan- 
dard deviations measured for the individual observations 
above, the constraints here from fitting the sets of obser- 
vations improved as little as a factor of 1.2 (PC mode 
photon-index) and as much as a factor of 3.5 (WT mode 
flux) . A greater improvement is achieved with WT mode 
than PC mode, which makes sense given that the indi- 
vidual WT mode observations have a large number of 
short (^ 2ks) observations. 

The probability of the data being constant can be 
estimated from the y 2 of each data set. For the PC 
mode observation sets, the reduced y 2 /^ for the flux is 
2.4/3 which corresponds to a 6.7% probability of being 
constant and for the power-law index the y 2 /i/ is 11/3 
(5.1 x 10~ 5 % probability). The sets of WT mode ob- 
servations have a y 2 /^ °f 2.3/8 (2.1% probability) for 
flux and a y 2 /;z of 6.1/8 (6.3 x 10 _6 % probability) for 
power-law index. So, in both modes, the fluxes are con- 
sistent with being constant (i.e. within 3<r), although the 
power-law indices are only consistent within a 5cr toler- 
ance. This could be due to unknown systematic sources 
of error or possibly low-level spectral variations, possibly 
due to the neglected blackbody component. Regardless, 
these variations are much lower than the ~ 30% pre- 
viously claimed (IRea et all 120051 ICampana et al.1 12007t 
IGotz et ahll2007lh 


3.2. Timing 

A phase-coherent timing solution for RXS J1708 has 
been maintained using Rossi X-ra y Timing Explorer 
since 1998; see iDib &; Kaspil (|2013l h In order to con- 
tinue to maintain a timing solution for RXS J1708, we 
began a monitoring campaign using the Swift /XRT on 
2011 July 28, overlapping with the RXTE campaign, un- 
til RXTE's demise in December, 2011. Monitoring ob- 
servations were typically 2-ks long. Barycentred events 
were used to derive a pulse time of arrival (TOA) for 
each observation. For a given observation, a TOA was 
obtained usin g a maximum like l ihood (ML ) method, as 
descri bed by ILivingstone et al.l (120091) and IScholz et al.1 
(|2012l! . The ML method compares a continuous model of 
the pulse profile, derived from taking aligned profiles of 
all the pre-glitch Swift /XRT observations, and creating a 
template composed of the first five Fourier components. 

These TOAs were fitted to a pulse arrival time model 


in which the phase, <j), at time t is given by: 

= 0o + v 0 {t - t 0 ) + -vo{t - t 0 ) 2 . (1) 

where <j> o, vq, and vq are the phase, frequency, and 
frequency derivative of the pulsar respectively at the 
reference epo ch tp. This was accomplished using the 
TEMP 02 (jHobbs et all [2006 1 pulsar timing software 
package. 

In Figure[5]we show the timing residuals for RXS J1708 
starting on 2011 July 28 and show the overlap between 
the RXTE and Swift monitoring epochs. The data are 
well fit by a single spin frequency and frequency deriva- 
tive (Tabic |2). However, we identified one notable timng 
event which we report as a new glitch. The event occured 
within 11 days (lcr uncertainty) of MJD 56019 with a de- 
caying Av/v = (8.3 ± 0.6) x 10- 7 . The glitch displayed 
an exponential recovery with a timescale of 111 ± 15 
days, 2.6 ± 0.3 times longer than the 43 ± 2 day decay 
ti me of t he othe r reported decaying glitch in the source 
()Dib et al.ll2008h . This glitch was accompanied by a Av 
of (1.4 ±0.3) x 10“ 15 Hz/s. 

4. DISCUSSION 

In this paper, we have reported on the flux and spec- 
tral properties of RXS J1708 over a ~ 10 year period 
from 2003 to 2013. We show that there is no significant 
flux variability and that only low-level spectral variations 
are seen. We have also presented an up-to-date timing 
solution and we report on a glitch that occured on MJD 
~ 56019. Below we compare our findings with previous 
results and discuss the significance of the lack of variabil- 
ity in RXS J1708. 

4.1. Flux variability of RXS J1708 

In this work, we do not use measured flux and spec- 
tral properties between different X-ray telescopes to con- 
strain the variability of RXS J1708. This is because 
cross-calibration between instruments onboard the Swift, 
Chandra, and X MM telescopes is such that the flux and 
spectral in dex can differ by up to 20% and 9%, respec- 
tively (e.g. iTsuiimoto et al.ll201lli . As seen in Figured! 
the X MM and Chandra observations are consistent with 
one another within those tolerances. Additionally, each 
Swift XRT mode (PC and WT) is considered separately, 
as the two modes are cross-calibrated only to within 10% 
in flux (A. Beardmore, private communication). 

Previous studies have claimed that RXS J1708 dis- 
played variability following glitches that occured between 

2002 and 2005. Us ing a multi-component blackbody plus 
power-law model, IGotz et all (120071) measure a low 1- 
lOkeV absorbed flux of ~ 3 x 10 -11 erg cm -2 s _1 for the 

2003 X MM observation and 2006 set of Swift observa- 
tions. They measured a higher flux for the 2004 Chandra 
observation and the 2005 set of Swift observations. Their 
highest flux measured is ~ 4.5 x 10 -11 ergcm~ 2 s -1 for 
the 2005 Swift set. This gives a total claimed variability 
of ~ 50%. 

Because we used a single-component model in Section 
IQ the values in Figure [T] are not directly compara- 
ble to those in IGotz et alJ ~ (j2007l ). For the Swift data, 
additional spectral components do not significantly im- 
prove the fit. However, additional components for Chan- 
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dra and X MM observations do provide a much bet- 
ter fit and so here we apply a blac kbody plus power- 
law model for direct comparison with IGotz et al.1 (120071) . 
With a multi-component model we measure 1-10 keV ab- 
sorbed fluxes of (3.83 ±0.04) x 10” 11 erg cm -2 s -1 for the 
XMM spectrum and (4.29 ± 0.08) x 10 -11 ergcm^ 2 s^ 1 
for the Chandra observation. This 11% discrepency in 
flux between the two obser vations is within the 20% 
cross-calibration error. As in IGotz et ail (120071 ). we find 
that the temperature of the blackbody component of the 
model is consistent between the two observations and 
is 0.46 ± 0.01 keV. For the photon index we measure 
2.63 ± 0.03 and 2.50 ± 0.07 for the XMM and Chandra 
observations, respectively. This com par es to T ~ 2.8 for 
both observations in IGotz et ahl (120071 ) . Reassuringly, 
the Chandra flux that we measure is higher than the 
XMM flux and the Chandra spectral i ndex is harder as 
found in cross-calibration studies (e.g. [Tsuiimot o et al.1 

iffl). 

In order to attempt to reproduce previous S wift re sults, 
for which only PC mode data were used (IGotz et al.1 
120071) . we processed the PC mode data from 2005-2007 
without using exposure maps and without removing or- 
bits with bad columns within 3 pixels of the center of 
the PSF. We found the same trend as in previous stud- 
ies: the flux of the 2005 set of observations was higher 
than the 2006 and 2007 sets and the flux of the 2007 
set was slightly higher than that of the 2006 set. We 
also observed an apparent correlation between the flux 
and power-law index. However, the level of variability in 
the three flux poi nts was only about 30% compared to 
~ 50% claimed in IGotz et ahl (120071 ) . Still, this is much 
higher than the < 10% that we find in our more detailed 
analysis. 

The lack of variability found here using soft X-ray 
imaging telescopes is consistent with what has been 
found by non-focusin g telescopes in other reg imes. Us- 
ing INTEGRAL data, den Hartog et al.l (|2008f ) found no 
significant variability in the hard X-ray flux or spectral 
in dex f or d ata sp anning from 2003 to 2006. With RXTE, 
iDib et all (j2008l) found that the pulsed count rate showed 
evidence for only low-level variability (< 15%) and they 
concluded that the glitches of RXS J1708 appeared to be 
“quiet”, i.e. unassociated with significant changes in the 
radiative properites of the magnetar. 

4.2. Radiative activity and Glitches in Magnetars 

Radiative activity in magnetars is almost always as- 
sociated with changes in timing beh avi or (e .g. glitches 
or increased timing noise: IDib fc Kas~pill2013D . Of the 26 
known magnetars and magnetar candidiates, only five 
have long-term (<10yr) phase-connected timing solu- 
tions that can be used to unambiguously detect glitches. 
These five magnetars are IE 1841—045, IE 2259±586, 
4U 0142±61, IE 1048.1-5937, and RXS J1708. Of 
the three glitches each that have been detected from 
IE 2259±586 and IE 1048.1—5937, five were radiatively 
loud, with the 2006 glitch of IE 2259±586 being the 
exception. The magnetars IE 1841—045, 4U 0142±61, 
and RXS J1708 have not displayed any significant 
flux increases associated with their glitches, although 
4U 0142±61 emitted short X-ra y bursts near the e poch 
of its 2006 candidate glitch (see IDib fc Kaspill2013l and 
references therein). 


It is therefore clear that glitches are not always ac- 
companied by radiative changes. Because changes in the 
magnetosphere would likely manifest as pulse profile or 
flux variations, it seems more likely that radiatively quiet 
magnetar glitches have their origin in the interior of the 
neutron star. If we assume that radiatively quiet and 
loud glitches have the same origin, a mechanism must 
exist to allow magnetars to exhibit prompt X-ray flux 
increases in some cases and no significant flux increases 
in others. 

One possible way to achieve both radiatively loud and 
quiet glitches in an interior model is to vary the depth at 
which the glitch-inducing event occurs. lEichler fc Chenil 
(119891 ) showed that if energy is injected into the crust of 
a neutron star it can travel outward, and manifest as a 
prompt outburst, or travel inward and heat the core of 
the neutron star. The direction of travel depends on the 
size and depth of the energy deposition. In the inward 
case, the heat is released slowly over a time scale of thou- 
sands of years. The flux decays of magnetars following 
prompt outbursts are indeed reasonably well modelled 
by crustal cooling (iLvuba rskv et all 120021: iScholz et al.1 
l2012t lAn et al.1 120131 : IScholz et a.1.1 l2013f ) . If the mecha- 
nism that causes glitches in magnetars injects energy at 
a shallow depth, a radiatively loud glitch would occur. 

An additional possible limit to the occurance of ra- 
diative outbursts from magnetars at glitch epochs is the 
predominance of neutr ino emission at high temperatures 
in neutron star crusts (Eichlcr & Cheng 1989; van Riper 
EMI). We expect neutron stars to have a limiting 
luminosity which occurs when the emission of neutri- 
nos dominates as a cooling mechanism over the emis- 
sion of photons. We would thus expect the brightest 
magnetars to be una ble to i ncrease t heir luminosity be- 
yond ~ 10 35 ergs -1 (iThompson fc Duncanl 11996 ). The 
five brightest magnetars, for which long-term timing so- 
lutions are available, have luminosities ~ 10 35 ergs _1 
(though see below for caveats on luminosity measure- 
ments). So, flux increases for these magnetars should 
either not occur or be small. Indeed, of the five, only 
IE 2259±586 and IE 1048.1—5937 h ave displayed signif- 
icant flux increases at glitch epochs (1 Woods et al.l 12004 
iGavriil fe Kaspil 12004 iTam et aD I2008D and those flux 
increases were much smaller than those from o utbursts 
observed in fainter transien t magnetars (e.g. Ilsrael et al.l 
l2007t IScholz fc Ka,snill20TTl) . 

However, magnetar luminosities are not well con- 
strained since the source distances are hard to deter- 
mine. There exist in the literature several disagree- 
ments in the distances to magnetars that lead to a dis- 
crepen cy of up to a fa ctor of ~ 30 in luminosity (e.g see 
lAn et all 120121 versus iDurant fc van Kerkwiikl 120061 for 
IE 1048.1—5937). Even when the distance is agreed upon 
there are discrepancies. For example, for RXS J17 08, the 
only distance estimation is from lDurant fc van Kerkwiikl 
(120061) but the 2-10 keV luminosity has been reported 
to b e as low as 4.2 x 10 34 ergs _1 (|Qlausen jkJ^aggjl 
2013) and as high as 1.4 x 10 35 ergs _1 dRea fc Esnositol 
201l|). From the model given by our best- fit mean flux 
and spectral indices, we get a 2-10 keV unabsorbed flux 
of 3.9 x 10~ n erg cm -2 s _1 , which corresponds to a 2- 
lOkeV luminosity of 6.8 x 10 3 4 ergs~ 1 , closer to tha t 
listed in the magnetar catalog (lOlausen fc Kaspil [2013). 


5 


Discrepancies such as this could be caused by the dif- 
ference in spectral models used or differences in the in- 
struments used to measure the flux (as mentioned above, 
X-ray detector cross-calibration can be discrepant up to 
20%). We therefore cannot say conclusively whether 
magnetar luminosity is inversely correlated with the size 
of radiative activity as d i scusse d here and as previously 
proposed in lPous fc Real (l2012h . 

If we do assume that fainter magnetars are able to 
have larger flux increases coincident with glitches, this 
suggests a rough luminosity order for the five brightest 
magnetars. RXS J1708, IE 1841-045, and 4U 0142+61 
have experienced only radiatively quiet glitches whereas 
IE 1048.1—5937 and IE 2259+586 have shown significant 
flux increases during some (or all for IE 1048.1—5937) 
of their glitches. That suggests that RXS J1708, 
IE 1841—045, and 4U 0142+61 are more luminous than 
the other two. 

Pulsars with higher B-fields are expected to be more 
luminous and have higher surface temperatures than pul- 
sars with lower magnetic fields because of energy depo- 
sition from the decay of their magnetic fields. Indeed, 
it has been shown that high-B radio pulsars are system- 
atically hotter t han similarily aged pulsars with lowe r 
magnetic fields dZhu et al.1 120111 : lOlausen et al.1 12013D . 
We may also expect that magnetar-like activity in such 
sources could arise due to energy from the magnetic field 
being deposited at shallow depths. Case in point, the 
high-B rotation-powered pulsar PSR J1846 — 0258 dis- 
playe d a magnetar-like outburst in 2006 dGavriil et al.l 
120081) . In recent years, two magnetars, SGR 0418+5729 
and Swift ,11822.3— 1606, were discovered with magnetic 
fields lower than several high-B rotation-powered pul- 
sars and have had clear X-ray outbursts (though it 
is unknown wh e ther or not they accom panied glitches; 
iRea et aD 120101 iLivingstone et al.l 1201 li b It is thus be- 
coming increasingly clear that high-1? rotation powered 
pulsars and magnetars are related and form a spectrum 
of objects rather than two distinct groups. Therefore, the 
mechanism that causes X-ray outbursts at glitch epochs 
could be active in all high-1? held pulsars. 

5. CONCLUSIONS 

We have presented an analysis of all of the Swift WT 
and PC mode data of RXS J1708 in the period 2005- 
2013. We show that the maximum variability for both 
the 1-10 keV X-ray flux and spectral index is constrained 
to < 10%. This is much less than claimed by previous 
studies and is consistent with the flux being constant. We 
also report on a newly discovered glitch at MJD~ 56019 
which has a fractional amplitude of Ai//zz = (8.3 ±0.6) x 
10 -7 , typical of magnetar glitches. 

The occurance of both radiatively quiet and loud 
glitches in magnetars, sometimes from the same source, 
shows that the mechanism that causes these glitches 
must be able to produce prompt flux increases in some 
cases and no signifcant increases in others. Here we 
have discussed the possibility that the glitches originate 
internally to the neutron star, with the deciding factor 
the depth of the energy deposition associated with the 
glitch. We note that these conclusions have been drawn 
from a sample of only five magnetars and therefore 
increasing the number of magnetars for which we can 
unambiguously detect glitches would be beneficial in 


answering the questions posed here. 
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Table 1 

Summary of Swift observations of RXS J1708 


Sequence 

Mode 

Observation date 

MJD 

Exposure time 

Set 

Set exp. time 

Set counts 




(TDB) 

(ks) 


(ks) 


00050701001 

PC 

2005-01-30 

53400.2 

2.3 




00050702001 

PC 

2005-02-02 

53403.0 

4.6 




00050702002 

PC 

2005-02-23 

53424.0 

2.0 

2005 

23.1 

10947 

00050701002 

PC 

2005-02-24 

53425.1 

11.9 




00050700006 

PC 

2005-03-23 

53452.2 

2.3 




00035318001 

PC 

2006-09-20 

53998.4 

2.7 

2006 

11.8 

5081 

00035318004 

PC 

2006-10-09 

54017.3 

9.2 

00035318005 

PC 

2007-02-25 

54156.3 

1.3 




00035318006 

PC 

2007-02-28 

54159.0 

1.8 




00035318007 

PC 

2007-03-05 

54164.9 

2.3 




00035318008 

PC 

2007-03-13 

54172.7 

1.2 

2007 

12.4 

5537 

00035318010 

PC 

2007-03-18 

54177.8 

2.0 




00035318012 

PC 

2007-03-23 

54182.6 

2.0 




00035318011 

PC 

2007-03-26 

54185.4 

1.7 




00035318013 

00090025001 

PC 

PC 

2008-02-23 

2008-05-13 

54519.2 

54599.0 

15.9 

5.0 

2008-PC 

20.9 

10169 

00090057001 

WT 

2008-04-02 

54558.0 

3.0 




00090057002 

WT 

2008-04-03 

54559.6 

2.1 




00090057003 

WT 

2008-04-04 

54560.5 

3.2 




00090057004 

WT 

2008-04-08 

54564.1 

1.2 

2008-1 

27.1 

26586 

00090057005 

WT 

2008-04-11 

54567.0 

3.8 




00090057006 

WT 

2008-06-05 

54622.2 

6.4 




00090057007 

WT 

2008-06-06 

54623.4 

7.5 




00090057008 

WT 

2008-08-13 

54691.3 

7.4 




00090057009 

WT 

2008-08-14 

54692.1 

1.4 

2008-2 

16.4 

15260 

00090057010 

WT 

2008-10-03 

54742.0 

7.0 

00090057011 

WT 

2008-10-10 

54749.7 

0.6 




00090057012 

WT 

2009-02-06 

54868.3 

2.2 




00090057013 

WT 

2009-02-08 

54870.1 

3.8 




00090057014 

WT 

2009-02-15 

54877.1 

12.7 

2009-1 

42.9 

42616 

00090057015 

WT 

2009-03-20 

54910.0 

16.5 




00090213001 

WT 

2009-04-26 

54947.1 

7.7 




00090213002 

WT 

2009-06-28 

55010.9 

8.7 




00090213004 

WT 

2009-09-02 

55076.2 

8.7 

2009-2 

22.3 

21748 

00090213005 

WT 

2009-10-11 

55115.0 

4.9 

00090213006 

WT 

2010-02-03 

55230.0 

8.6 




00090213007 

WT 

2010-02-04 

55231.7 

5.4 

2010 

23.8 

23038 

00090213008 

WT 

2010-03-25 

55280.0 

9.8 




00035318014 

WT 

2011-07-28 

55770.3 

0.9 




00035318015 

WT 

2011-08-04 

55777.4 

1.0 




00035318016 

WT 

2011-08-11 

55784.0 

2.3 




00035318017 

WT 

2011-08-18 

55791.2 

1.9 




00035318018 

WT 

2011-08-25 

55798.3 

2.0 




00035318019 

WT 

2011-09-01 

55805.4 

2.1 




00035318020 

WT 

2011-09-08 

55812.5 

2.2 




00035318021 

WT 

2011-09-15 

55819.2 

2.2 

2011 

27.8 

27744 

00035318022 

WT 

2011-09-22 

55826.1 

2.3 




00035318023 

WT 

2011-09-29 

55833.9 

2.4 
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Table 1 — Continued 


Sequence 

Mode 

Observation date 

MJD 

Exposure time 

Set 

Set exp. time 

Set counts 




(TDB) 

(ks) 


(ks) 


00035318024 

WT 

2011-10-06 

55840.8 

2.3 




00035318025 

WT 

2011-10-13 

55847.0 

1.6 




00035318026 

WT 

2011-10-20 

55854.3 

0.9 




00035318027 

WT 

2011-10-22 

55856.3 

1.8 




00035318028 

WT 

2011-10-27 

55861.1 

2.0 




00035318029 

WT 

2012-01-25 

55951.5 

2.0 




00035318030 

WT 

2012-02-01 

55958.3 

2.1 




00035318031 

WT 

2012-02-08 

55965.8 

2.2 




00035318032 

WT 

2012-02-15 

55972.2 

0.3 




00035318033 

WT 

2012-02-22 

55979.8 

0.4 




00035318034 

WT 

2012-02-29 

55986.4 

2.2 




00035318035 

00035318036 

WT 

WT 

2012-03-07 

2012-03-16 

55993.3 

56002.6 

2.2 

2.2 

2012-1 

21.0 

21056 

00035318037 

WT 

2012-03-21 

56007.5 

1.0 




00035318038 

WT 

2012-03-28 

56014.9 

1.6 




00035318039 

WT 

2012-04-13 

56030.2 

1.6 




00035318040 

WT 

2012-04-26 

56043.1 

1.5 




00035318041 

WT 

2012-05-10 

56057.3 

1.3 




00035318042 

WT 

2012-05-25 

56072.0 

0.5 




00035318043 

WT 

2012-06-07 

56085.6 

2.0 




00035318044 

WT 

2012-06-22 

56100.2 

1.5 




00035318045 

WT 

2012-07-05 

56113.4 

0.6 




00035318047 

WT 

2012-07-15 

56123.5 

1.1 




00035318048 

WT 

2012-08-16 

56155.0 

1.7 

2012-2 

10.8 

11423 

00035318049 

WT 

2012-09-06 

56176.1 

2.0 




00035318050 

WT 

2012-09-27 

56197.1 

1.7 




00035318052 

WT 

2012-10-23 

56223.4 

0.3 




00035318053 

WT 

2013-01-24 

56316.3 

1.0 




00035318054 

WT 

2013-02-13 

56336.8 

1.9 




00035318055 

WT 

2013-03-06 

56357.1 

2.0 

2013 

7.1 

7380 

00035318056 

WT 

2013-03-27 

56378.6 

0.7 




00035318057 

WT 

2013-03-31 

56382.2 

1.5 





Table 2 

Timing Parameters for RXS J1708. 



Parameter 

Value 

Observation Dates 


28 July 2011 - 29 May 2013 

Dates (MJD) 


55770.396 - 56441.770 

Epoch (MJD) 


56000.000 

Number of TOAs 


61 

* (O 


0.090851264(3) 

v (s 2 ) 


-1.638(3) x 10 -13 

Glitch 

Glitch Epoch (MJD) 


56019(11) 

A v d (s- 1 ) 


7.5(5) x 10" 8 

T d (days) 


111(15) 

A v (s -2 ) 


1.4(3) x 10” 15 

RMS residuals (ms) 


229.07 

x 2 /*' 


64.94/55 

Numbers in parentheses 

are TEMP02 reported 1 cr uncertainties. 




53000 53500 54000 54500 55000 55500 56000 

MJD 


Figure 1. Top Panel: Absorbed 1-10 keV flux of RXS J1708 over a ~10-yr period. Note that the zero on the y-axis is suppressed. 
Bottom Panel: Photon indices from fitting a power-law model to the 1-10 keV spectrum. Grey points are from spectral fits to individual 
observations. Black triangles are sets of Swift PC mode observations, and white circles are sets of Swift WT mode observations (see Table 
[I] for definitions of sets). XMM and Chandra observations are labelled. The dark grey bands represent the 90% error in the mean and 
the l ight grey bands represent the level of previously claimed variabilty (~ 50% in 1-10 keV flux and ~ 30% in spectral index; IGotz et al.l 
120071) . The solid vertical lines represent the epochs of glitches and the dashed lines indicate the epochs of glitch candidates. All error bars 
are 90% confidence intervals. 
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on 

cn 

73 

P 

• i-H 

on 

<n 

pc; 


on 


on 

7d 

• i— I 

on 

pc; 



55700 55800 55900 56000 56100 56200 56300 56400 56500 


Modified Julian Date 


Figure 2. Timing residuals, the difference between the predicted and measured TOAs for the timing model shown in Table [2] The top 
panel shows residuals before fitting for a glitch, and the bottom panel after. In both panels, open circles indicate data from RXTE , and 
black triangles indicate Swift. The vertical dashed line represents the glitch epoch and the grey band represents the uncertainty in that 
epoch. 


